DMRG analysis of the SDW-CDW crossover region 
in the ID half-filled Hubbard-Holstein model 
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■ Abstract. In order to clarify the physics of the crossover from a spin-density- wave (SDW) 

Mott insulator to a charge-density-wave (CDW) Peierls insulator in one-dimensional (ID) 
systems, we investigate the Hubbard-Holstein Hamiltonian at half filling within a density 
• matrix renormalisation group (DMRG) approach. Determining the spin and charge correlation 

exponents, the momentum distribution function, and various excitation gaps, we confirm that 
an intervening metallic phase expands the SDW-CDW transition in the weak-coupling regime. 
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The Hubbard-Holstein model (HHM) pQ is archetypal for exploring the complex interplay 

of electron-electron and electron-phonon interactions especially in quasi- ID materials, such as 

halogen-bridged transition metal complexes, charge transfer salts, or organic superconductors [2]. 

It accounts for a tight-binding electron band (ex 2t), an intra-site Coulomb repulsion between 

electrons of opposite spin (ccu = U /At), a local coupling of the charge carriers to optical phonons 

(oc A = g 2 uo/2t), and the energy of the phonon subsystem in harmonic approximation (oc oJo/t): 
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Here c ia (c ia ) creates (annihilates) a spin-u electron at Wannier site i of an ID lattice with 
N sites, riic = c ia c ia , and b\ (6J are the corresponding creation (annihilation) operators for a 
. dispersionless phonon. We consider the case J2 icr ni a = 1 hereafter, and take t as energy unit. 

Based on exact diagonalisation data for the staggered static spin/charge structure factor, 
Sa/p{l) = jflt2j£ eiq ^(( n j't ^ n jl)( n £t i ^ has been argued that the HHM shows a 

^ . crossover between Mott and Peierls insulating phases near u/X ~ 1 |3j. But this only holds in 

the strong-coupling adiabatic-to-intermediate phonon frequency regime. Later on the ground- 
state phase diagram of the HHM was explored in more detail, also for weak interaction strengths 
and large phonon frequencies. In this regime, variational displacement Lang-Firsov [1], stochastic 
series expansion QMC [5], and DMRG [6] methods give strong evidence that, if A is enhanced at 
fixed u and uq, the SDW-CDW transition splits into two subsequent SDW- metal and metal-CDW 
transitions at A c i and A C 2, respectively (see fig. [H dashed and dot-dashed lines). Very recent 
DMRG data indicated that in the anti-adiabatic regime of very large phonon frequencies the 
metallic phase might be even more extended than the one obtained by QMC and is subdivided 
into regions with a normal ID metallic (I) and a bipolaronic-liquid (II) behaviour [7]. 

In this work, we will re-examine the weak-coupling SDW-CDW transition regime by 
calculating the ground-state properties of the HHM in the framework of a large-scale numerical 
(boson pseudo-site) DMRG approach supplemented by a finite-size scaling analysis [8] . 



Figure 1. Qualitative phase diagram of the ID Hubbard- 
Holstein model. Given that in the half-filled HHM model 
the ground-state is metallic at u = for ojq > provided 
that g < g c , it was proposed that this metallic phase 
continues to exist between the SDW and CDW states for 
u > [H El El E]- With increasing ujq the region of the 
intervening metallic state increases, and the tricritical point 
u m moves to larger u [5|- The SDW state shows no long- 
range order (LRO) and is characterised by a vanishing spin 
gap A s but a finite charge gap A C1 , whereas the CDW phase 
exhibits true LRO and A s = A C1 > 0. 



Characterising the SDW-CDW-intervening metallic phase of the HHM, we presume that 
for the metallic state a Tomonaga-Luttinger-liquid (TLL) description holds. In the TLL 
picture, nonuniversal coefficients, K p and K a , determine the decay of correlation functions 
and therefore can be used to identify the properties of the TLL phase [9], but also the phase 
boundaries to the insulating states [S]. In practice, we can extract the TLL correlation exponents 
from the slope of the corresponding structure factors in the long- wavelength limit HUj : 
K P /a = vrlim^o S p/r7 (q)/q, where q = 2ir/N for N ->■ oo. 

Specifically, K p > 1 [K p < 1) corresponds to attractive (repulsive) charge correlations in 
the TLL and K p = signals an insulating phase. Hence K p jumps from 1 —?■ at the metal- 
SDW/CDW transitions. The spin exponent takes the value K a = in a spin-gapped phase and 
K a = 1 everywhere else in the thermodynamic limit [11] . For finite systems the situation is more 
involved, in particular for the spin exponent K a . First, the convergence K a — > is slow-going as 
N — > oo in the spin-gapped phase. Second, logarithmic corrections prevent K a — > 1 in the spin- 
gapless (SDW) phase. On the other hand, these logarithmic corrections vanish at the critical 
point, where the spin gap opens, and we can utilise that K a (K p ) crosses 1 from above (below) 
at some A Cl (as the electron-phonon coupling increases for fixed u), in order to determine the 
SDW-metal phase boundary itself. Increasing A further, K p should cross 1 once again, this time 
from above, at another critical coupling strength, A C2 , which pins the metal-CDW transition 
point down. 

Figure [2] corroborates this scenario for the anti-adiabatic regime of the HHM. The two critical 
values A Cl and A C2 are in accord with the phase diagram obtained by QMC [5]. K a < 1 and 
K p > 1 earmark the intervening metallic phase. In terms of the TLL framework, a metallic phase 
with K p > 1 exhibits dominant superconducting correlations. Recent DMRG calculations of the 
the (s-, p-, and cf-wave) superconducting correlation functions of the half-filled HHM indicate, 
however, that these correlations are only sub-dominant against CDW correlations [6], while QMC 
investigations attributed the K p > 1 to finite-size effects and suggest that K p (N — > oo) = 1, 
i.e., superconducting and CDW correlations are exactly degenerate. 

Here we inspect the finite-size scaling of the spin and single-particle charge excitation gaps, 
A s (A0 = E (l) - E (0) and A Cl (N) = E+(l/2) + £^(-1/2) - 2E (0), respectively, as well as 

that of the two-particle binding energy A b (N) = E%~(Q)+Eq(0)-2E%(-1/2), where E { Q L±) (S z ) 
is the ground-state energy at or away from half-filling with N e = N ± L particles in the sector 
with total spin-z component S z . The left panel of fig. [2] shows that both spin and charge 
gaps open at A Cl (but there is no LRO). For u < u m , the transition at A Cl seems to be of 
Kosterlitz-Thouless type, i.e. just above A Cl the gaps are exponentially small and therefore their 
magnitude is difficult to determine. In this region, denoted by (I) in fig. [H we find A C1 ~ A s , 
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Figure 2. DMRG results for the half-filled HHM with u = 0.25 and ojq = 5. Left panel: charge 
(K p , filled circles) and spin (K a , open squares) TLL exponents as functions of A [left-hand axis 
of ordinate]. Open triangles up (triangles down) give the spin (charge) gap A s (A Cl ) [right- 
hand axis of ordinate]. Data shown are extrapolated values for the infinite systems, using open 
boundary conditions. Right panel: Finite-size scaling of the binding energy Af, for different A; 
lines are polynomial fits. In the numerical calculations we use up to five pseudo-sites and keep 
2400 density-matrix eigenstates; then, for all parameters studied, the local boson density is less 
than 10 -8 and the discarded weight is smaller than 10 -9 . 



and the binding energy Af, is also extremely small, or maybe even zero (see triangles up, right- 
hand panel of fig. [2]). As A increases, we obtain a (smooth) crossover to a metallic regime with 
a noticeable two-particle binding Af, < (region (II) in fig. [1]), where A Cl ~ A s . This is in 
accord with the very recent findings of Ref. [7|, where a subdivision of the metallic phase into a 
weakly renormalised TLL (I) and a bipolaronic liquid^ (II) was suggested. In the latter phase, 
the two-particle excitation gap A C2 (iV) = Eq + (0) + Eq~(0) — 2Eq(0) was shown to scale to zero. 
In the CDW phase, which typifies a bipolaronic superlattice at large phonon frequencies, we 
have, besides A s = A C1 > 0, A C2 > and Af, < 0, whereas in the SDW state A C2 > but 
Af,(iV — > oo) — > 0. While the basic scenario discussed so far persists in the adiabatic regime, 
the metallic region shrinks as the phonon frequency ojq becomes smaller [5J [7]. Furthermore, 
the CDW state rather behaves as a normal Peierls insulator and consequently there is a weaker 
tendency towards bipolaron formation in the metallic state for small A, and u < u m . 

Finally, let us investigate the behaviour of the momentum distribution function, n a (k) = 

jj X^z=i cos (Mi — 0) ( c ] a c i o-)> ^ = 27rm/./V, m = 0, . . . , N/2, which can be obtained by DMRG 
for a system with periodic boundary conditions. 

Figure [3] shows the variation of n(k) for weak (circles) and intermediate (stars) Hubbard 
interactions in the SDW (a), TLL (b)-(c), and CDW (d) phases. The momentum distribution is 
a monotonously decreasing function as k changes from the centre (k = 0) to the boundary 
of the Brillouin zone (k = ir). Since we consider the weak-coupling regime, n(k) is only 
weakly renormalised away from the Fermi momentum kp. For a ID TLL, instead of the 
Fermi liquid typical jump of n(k) at kp, one finds an essential power-law singularity [9], 
corresponding to a vanishing quasiparticle weight Z = 0. For finite TLL systems, the difference 
A = n(/cF — 5) — n(kp + 5) is finite (with 5 = tt/N = tt/66 in our case), and rapidly decreases 
with increasing couplings A, u. Approaching the insulating SDW/CDW states n(k) becomes a 
smooth curve, i.e. the singularity vanishes and A — > 0. At very large A, the system develops a 



1 Note that the polaronic two-particle bound states are not necessarily small (i.e. on-site). 



Figure 3. Momentum distribution of the 
half-filled HHM in the anti-adiabatic regime 
(cjo = 5). Open (closed) circles give DMRG 
results at u = 0.25 (u = 1) for a system 
with N = 66 sites and periodic boundary 
conditions. The occupation of fermionic 
states carrying momentum k is given by 
n(k) = \ ^2 a n a {k), and we have kp = ir/2, 
rik F = 1/2 for the half-filled band case. In 
the intermediate metallic phase, n(k) exhibits 
a power-law singularity at Hf [see panels (b) 
and (c)]. At weak and strong electron-phonon 
couplings insulating SDW [panel (a)] and 
CDW [panel (d)] are realised, respectively. 



"perfect" CDW with n(k) = 1/2 for all momenta k. 

To summarise, we validated the existence of an intervening metallic phase in the SDW-CDW 
transition regime of the ID half-filled Hubbard-Holstein model for u < u m by large-scale DMRG 
calculations. Spin and charge gaps open exponentially slowly at the SDW-TLL transition point, 
A Cl , but no long-range order develops. A Cl (u,u;o) can be determined from the "l"-crossing of 
the spin and charge TLL parameters. In the TLL, the momentum distribution function exhibits 
a power-law singularity at kp. As the electron-phonon coupling increases, a crossover to a 
bipolaronic metal, indicated by negative binding energy, takes place, before the systems enters 
the long-range ordered insulating CDW phase at a second critical coupling A C2 . We would like 
to point out that fixing the metal-CDW phase boundary quantitatively is a difficult issue. 
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